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We display several examples of how fields with different limiting velocities (the "speed of light") 
at a high energy scale can nevertheless have a common limiting velocity at low energies due to the 
effects of interactions. We evaluate the interplay of the velocities through the self-energy diagrams 
and use the renormalization group to evolve the system to low energy. The differences normally 
vanish only logarithmically, so that an exponentially large energy trajectory is required in order to 
satisfy experimental constraints. However, we also display a model in which the running is power- 
law, which could be more phenomenologically useful. The largest velocity difference should be in 
system with the weakest interaction, which suggests that the study of the speed of gravitational 
waves would be the most stringent test of this phenomenon. 
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1. INTRODUCTION 

Many physical systems yield wave-like solutions which 
satisfy the wave equation with a speed of propagation Ci, 



9t2 



(x,t) =0 



(1) 



which is also the massless Klein-Gordon equation. To 
leading order, the Lagrangian of any such field obeys a 
Lorentz-like symmetry of Lorentz transformations scaled 
with the limiting speed c^, even if the underlying sys- 
tem does not have that invariance. However, if there are 
multiple fields, they will in general have different lim- 
iting velocities, and there will not be a global Lorentz 
symmetry. If the fundamental interactions are emergent 
phenomenon Jrom an underlying theory without Lorentz 
invariance [IhSj], we might expect that particles would dis- 
play different limiting speeds. 

In this paper we show how interactions between the 
fields can lead to a universal limiting velocity, i.e. the 
speed of hght, at low energies. We calculate how the 
different fields influence each other's propagation veloc- 
ity through the self-energy diagrams, and then use the 
renormalization group to scale the results to low energy. 
Using several examples we show that the condition of 
equal velocities is the low energy endpoint of renormal- 
ization group evolution^ . A heuristic explanation for this 
is that because fields can split into other types of fields, 
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^ During the course of this work we found that this general ap- 
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ization group running that we describe is also related to the study 
of the renormalization of Lorentz- violating electrodynamics @| . 



the propagation velocity of one field approaches that of 
the related other fields. 

This result could be useful if the fields of the Standard 
Model are emergent from an underlying theory that is 
not Lorentz invariant. Of course, Lorentz invariance is 
conventionally taken as one of the foundational princi- 
ples underlying all our fundamental interactions. How- 
ever, the Weinberg- Witten Q theorem is usually inter- 
preted as telling us that non-Abelian gauge bosons and 
gravitons cannot be emergent fields arising from any un- 
derlying Lorentz invariant non-gauge theory. All known 
examples [T| satisfy this property. Therefore if the idea of 
emergent fields has any application in the fundamental 
interactions it appears required that Lorentz invariance 
is also emergent. 

Our results show that a universal limiting velocity can 
be an emergent property in the low energy limit. How- 
ever, in general the difference in velocity runs towards 
zero only logarithmically. This means that the underly- 
ing scale of emergence needs to be exponentially far away, 
making it difficult to test any feature of that theory which 
is power-suppressed. For example, we estimate that the 
scale where differences in the velocity are of order 10% 
would be beyond 10^° GeV. Because of this feature we 
propose a model that produces much faster power-law 
running. The model involves a large number of fields 
which accelerate the running. Another consequence of 
the running speeds is that observable differences in the 
velocities would be greatest if the interactions are the 
weakest. This suggests that the measurement of the ve- 
locity of gravitational waves would be the most sensitive 
test of this aspect of emergence. 

This paper has the following structure. In the next 
section we give some general comments on our proce- 
dure. Then in Sec. 3, 4, 5 and 6 we calculate the beta 
functions for Yukawa theories, electrodynamics, Yang- 
Mills and mixed theories respectively. All cases yield 
beta functions such that the limiting velocities run to- 
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wards each other at low energy. In Sec 7, we analyze the 
general effects of logarithmic running and address the 
phenomenological constraints. Because of the difficulties 
posed by logarithmic running, we display a model with 
power-law running in Sec. 8. We close with a summary 
and discussion. Some of the more technical details are 
described in a pair of appendices. 



rections, while those of the couplings require vertex cor- 
rections as well (see Fig. [T]). 




2. SETUP 

We assume that different species of fcrmions, scalars 
and gauge fields emerge at some UV scale with different 
speeds of light. In condensed matter systems, phonons 
and magnons do not propagate at the same speed. Sim- 
ilarly, the same behavior is expected to carry on in an 
emergent theory of nature. In the absence of any form of 
interactions between particles, their speeds are expected 
to be frozen as we run down to the IR. However, these 
particles are interacting due to Yukawa and gauge forces. 
Hence, the total Lagrangian of such a system will be 
given by the sum of kinetic and interaction terms with 
certain bare coefficients specified initially at the UV. The 
parameter space of the system is spanned by the dif- 
ferent speeds and interaction strengths. According to 
the principle of self-similarity and Wilsonian renormal- 
ization, the same Lagrangian will continue to describe 
the system at different energy scales provided that we 
replace the bare parameters with the renormalized ones. 
This can be achieved by integrating out the high momen- 
tum modes as we run down from UV to IR. Quantum 
loops are sensitive to high momenta, and hence can be 
used to track the evolution of trajectories of the different 
speeds and interaction strengths in the parameter space. 
The evolution of these trajectories are encoded in the f3 
functions that are given by the Gell-Mann Low equations 



Pigt) = = f{9j} : 



where is the mass scale we introduce in dimensional 
regularization. 

In theories with a universal limiting velocity, the 
Lorentz symmetry prevents the renormalization of the 
speed of light, and one can set c = 1 as a definition of 
natural units. However, if different species carry differ- 
ent limiting velocities, Ci, then these parameters also get 
renormalized and must be treated in the same manner 
as coupling constants. They carry a scale dependence 
through the renormalization procedure, and also gener- 
ate their own beta function. We exploit this property to 
study the running of the limiting velocities. 

Throughout the paper we use dimensional regulariza- 
tion (dim-reg). The high energy part of the quantum 
loops can be isolated by retaining only the 1/e pieces 
that arise upon using dim-reg. Finally, we notice that 
our treatment is limited to one-loop corrections, and that 
the /? functions of the speeds require only self-energy cor- 



FIG. 1: The self-energy and vertex diagrams. Only self- 
energies will contribute to the running of the speeds, while 
the vertex is needed for the running of the coupling strength. 



3. YUKAWA INTERACTIONS 

We consider a two-species system, namely scalars and 
fermions, having different speeds of light at the UV and 
coupled through Yukawa interaction. The Lagrangian 
density reads 



t>r ■ d(j)r — g(t>r4'r^r ■ 



(3) 



where the subscript r denotes the renormalized values of 
the fields. The momentum-space propagators for scalars 
and fermions are given by 



Sf{p\p) 



(p0)2 _ ^2p2 

i 



(4) 



(2) The self-energies of fermions and scalars are respectively 



^Dk{p° -q\p-q) , (5) 



and 



iH(p",p) 



x5/(/ + g°,g+p)] . 



(6) 



In the following, we will be interested only in the di- 
vergent pieces of ([5]) and ([6]) . The integral ([6|) is trivial to 
perform upon using the substitutions h!^ — / Cf, k = q, 
po — pO ^^^^ gj-^^ p — Then, one readily finds 



.0\2 



finitej . (7) 
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On the other hand, the integral ([5]) is more involved and 
it needs a bit more attention. Using the substitution 
— q^/cf, and k = qwc find 



f d'^k 1/ 



{p° /cb- Cfk° /cb) - [k ~ p 



7"/" - 7 • / 



(8) 



where the integrals I*^ and / are given by (the details are 
in Appendix A) 



2cb 



(47r)^c2(l + a)2 



— h finite | 



and 



ip 2a(l + 2a) f2 
(inf 3(l + a)2 



finite I 



(9) 



(10) 



where a = Cb/cf. 

Now, we move to the vertex correction which, to one- 
loop order, reads 



zgG = {~^gf [ ^sM-q',p,-q) 
J (27r) 

xSf{p°-q°,Pi^q)Db{q°,q). 



(11) 



Using the change of variables q'^/cf = K'^, q — K, 

Pl/cf = P2, P2 = P2, P? = P?/cf, and pi = Pi, and 
retaining only the divergent part of the integral we ob- 
tain (see Apendix A) 



zgG- 



(47r) cjcb{l + a) 



— h finite I 



(12) 



At this point, the total Lagrangian including the one- 
loop effect is 



Cq — Cr + Cc 



(4^)'(l + a)2 
a(l + 2a) 



CbCf 



3cl 
9' 
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9o4'rdQ(f'r — 94>r ' d4>r 
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— 2 1 + a \ g 



(4^)' Cbc} 

where is the counter Lagrangian 



(13) 



'-do(l)rdQ(t>r 



'-b'-'fr 



(14) 



At this point, we can read off the different 5s that are 

required to absorb the infinities. Furthermore, we define 

1/2 1/2 
the bare fields 0o — <pr, and ipa = 'i'r, bare 

speeds Cf„, and Cf,^ and bare coupling go such that the 

Lagrangian density reads 

1 



(15) 



Comparing and (ITSt we find 



CbZ' 



-1/2 ryl/2 



ry ry-l/2ry-l e/2 

50 = 9ZgZ^ ' Z^ ^i' 



(16) 



where Z = I + 6 

= 1 — 



29' 



-4> 



1 



(47r) C6(c/ -|-C6)2 V e 
2.92 /2 



^/ = 1 
Zb = 1 



= 1 



(4^)' 4 Ve 

2g2(c^. + 2cb) 
3(47r)^c/Cb(c/ -f C6)2 
2g2 /2 



2g' 

(47r)^ CfCb{cf + Cb) 



(17) 



To proceed, we regard all the renormalized quantities 
above as functions of the scale /Lt that occurs in dim-reg. 
Then, we differentiate the system in Eq. [T6]w.r.t /i and 
solve simultaneously for l3{g), f3{cb) and l3{cf) to find 



P{cb) = 



9' 



[icbc) + 2clcf + 4 + 



Sn^Cbcjicf + CbY 



9' {^l - 4) 

^TT-^Cbc) ' 

,g2(c/ - Cb) 
6T:^Cb{cf + CbY 



(18) 



Notice that the /? functions of Cf, and c/ do not depend 
on the vertex correction Zg. Finally, we calculate the /3 
function of the ratio a ~ Cb/cf to find 



o, N /3(cb) Cb , , 
p{a) = — -^P{cf) 
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487r2 



Cf,c2(l -|-a)2 



(19) 



from which we see that Cf, = c/ is an IR attractive 
line. We can also see that by studying the Jacobian 
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J = dp(ci) / dcj, for i, j = c/, Cb at the fixed line c/ = Cb- 
The eigenvalues of J are {0, 7g^/247r^ci^}. The positiv- 
ity of the second value ensures that c;, = c/ is an IR 
attractive fixed line. 

We have seen the existence of an attractive IR fixed 
line corresponding to a common limiting speed. We will 
address more details about the running in Sec. 7. 



4. NON-COVARIANT ELECTRODYNAMICS 

In this section, we study the RG fiow of the limiting 
speeds of fermions and photons. The non-covariant and 
gauge invariant Lagrangian density reads 

1 



The counter Lagrangian reads 

= ^cgaugc+iSz^i>rdo'y°'^r--i3zfCf-ipr'l-d-lpr, 

(25) 

and £c gauge is the counter term for the gauge sector. 
Then, from (l^5t and (P^ . and after using the properties 
of 7 matrices, wc can immediately read and Z f 



Z^ — 1 — 



(20) 



{Any c,{cf + 
^ _ ^ 2eHcl+c})i2cg+Cf) /2\ 
^ Si4nfcfCgicg + Cf)^ [ej - 

Now we come to the counter terms in the gauge sector. 
A general counter term written in the momentum space 
takes the form 



where Frf,u = - duArf,, and 9^ = {do,Cgd), and 

Cg is the photon speed. The photon propagator in the 
Feynman gauge is given by 



(21) 



To find the photon and fermion self-energies, it 
proves easier to write the interaction Lagrangian in 
the form £/ — — eCpy-0rA(f7'^V'r, where c^v = 
diag(cg, — c/, — c/, — c/). Hence, the fermion self-energy 
is 



7"/" - 7 • / 



(22) 



where /q and I\ are given in ([9]) and (jlOp after replacing 
Cf, with Cg. While the photon self-energy is given by 



+S'^ {Sa {p'Y 

-2AirSACgP°p'^Aor ■ 



A 



1 r 

(26) 



One can show that all the infinities in can be ab- 
sorbed using 6a and Sg^ 



3(47r)'c/ \e 



Zgs — 1 - 



4 e^cf (2 



3(47r) ^g 



(27) 



where as usual Z — 1 + 6. We write £c gauge (p) in the 
compact form 



c gauge 



{p)^A^rM^'-A^, 



(28) 



with 



— — jtr [YSf {q°,q) 
(27r) 

7^5/(p" + g°,j?+-f)] . (23) 

Using the substitution fcp — q'^/cf, k — q, Pq — p'^/cf, 
and P — p we can put IIq,^ in a standard integral form. 
Hence, 



i^ap {p",p) 



Ai 



6 Cajy^Pn 



3 (47r)" c/ 
'2' 



(24) 



where P = [p'^ /cf,p). Explicit calculations shows that 
"P^IIq^ = 0, where V"' = (jP.Cgp), and hence Hq,^ is 
gauge invariant as expected. 



-Cg6Ap"p^ , 



(29) 



It is trivial to see that VaM"^ = 0, and hence Map is 
gauge invariant. 

Now, defining the bare fields ipo — Z^^^Tpr, Aq — 

Z\Zg'g^'^A^, and Aq — Z^^A^, and bare speeds c/o, 
and Cgo the Lagrangian density reads 



Co = -^-^OtiiyFl^" + iipodoj^^o - ic/oV'o7 ' ^V'o ■ 



(30) 
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The bare gauge field Lagrangian in the momentum space 



is given by Aq^MI^^Aq^, and 



(31) 



The relations between the bare and renormalized speeds 
are 



CgO = CgZ^ ZJ^ , 

from which we obtain 



(32) 



^^'''^ " 3(4.)^ c,c, ' 

(c/ - cg) + 3c/c<, + c2) 

/?(c/) = ^77-2 V— ^2 ^-(33) 

3 (47r) Cg (C/ + Cg) 

These /? functions have the same structure as in the case 
of Yukawa interactions, and we immediately conclude 
that = Cg is an IR attractive line. 



5. NON-COVARIANT YANG-MILLS THEORIES 

In this section we generalize the results of QED to the 
case of non-abelian gauge theories. We take the gauge 
group to be SU{N), and the fermions in the fundamental 
representation 

— free 

+gc^,AA7YAt'' - gCgr'-d,A-,^AfA^- 

(34) 



2 2 J'eab ^ecd Aa Ab ack \d\ 



where g is the coupling constant, t° are the group gener- 
ators and /"^"^ are the group structure constants. £rfroe 
is the free part of the Lagrangian 
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—icfijjrd-^ipr, (35) 
where as in the case of QED — dfj^A^:^ — d^^A'^^, 

df_L = {do,Cgd), and Cg is the gauge boson speed. 

The fermion self-energy is identical to the case of QED, 
one just includes the quadratic Casimir operator in the 
fundamental representation C2{N) — {N'^ — 1)/2N into 
Eq. [mto find 



^ 2C2iN)g'i3c}-cl) /2\ 
(47r)2 cg(c^ + cg)2 [ej ' 



Zf 



3(47r)'c/Cg(cg + c/)2 \e) 



In calculating the gauge boson self-energy 11°^ one en- 
counters, in addition to the fermion loop, gauge boson 
and ghost loops 



. 5C2(G)g^ ^ ^2 N 



3 (47r)' 
AC{N)g^ 



3cf 



srab 



(36) 



where C2(G) = N, and C{N) — 1/2 are group factors, 
7"^ = (^p^jCgp) and — {p^/cf^p). As we did in 
QED, the infinites can be absorbed into the counter term 
A°^M'^''yl°,, where M^"' are given in Eq. [29l Hence, we 
find 



1 



ZgB — 1 



AC{N)g^ , 5C2(G).g2 



3(47r)^c/ 3(47r)^Cg 



^C{N)g'cf , 5C2(G)g 



3 {AttY cl 3 {AttY c 




(37) 



Gluon loops will not modify their own propagation speed, 
due to the Lorentz-like symmetry of that sector when 
considered in isolation. This is visible in the formulas 
above. Since /3(cg) cx {Zg^ — Z^), the gauge bosons and 
ghosts contributions cancel in obtaining l3{cg). Overall, 
the /3 functions read 



PiCg 



AC{N)g' 



3(47r) 



8C2iN)g^ (c/ - Cg) (44 + ScfCg + cfj 



3(4^)' 



Cg (Cf 



(38) 



which, apart from group factors, are identical to the QED 
case. 



6. EMERGENCE OF LORENTZ SYMMETRY IN 
A MIXED SYSTEM 

The emergence of a universal Lorentz Symmetry in the 
above examples is intriguing to explore a more general 
setup consisting of multi-species and/or mixing between 
fermions, bosons and gauge fields. Before delving into 
the most general case we derive a general formula that 
enables us to calculate the f3 functions of such complex 
systems. This is done in Appendix B. 



6.1. Yukawa-Electrodynamics 

Now, let us consider the more general case of Yukawa- 
electrodynamics. In this theory a fermion couples to a 
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scalar through Yukawa interaction, and minimally to a 
U{1) gauge field. The scalar is neutral under the [/(I) 
field. We assume that the fermion, scalar and gauge field 
all have different speeds of light, c/, Cb and Cg respec- 
tively. This is the simplest generalization of the above 
cases. The scalar and gauge field self-energies are iden- 
tical to their expressions in Yukawa and QED sections, 
while the fermion self-energy is the sum of the contribu- 
tions from the scalar and gauge field. The calculations of 
the corresponding Z renormalizations are very straight 
forward, and can be obtained directly from the previous 
two sections. Thus, Z^, Z^, Za, and are given by 
their expressions in Eq. [TTland Eq. [27lresDectivelv. while 



Z,i, 



2e^(3c^-cg) 



Zf = 1 



2e\cl + c}){2cg + Cf) 
3 (47r)^ c/Cg(cg + c/)2 

2g^{cf + 2cb) \ (2 



3 (47r) CfCb{cf + Cb)2 



(39) 



The relations between the bare and renormalized quan- 
tities are given as before 



c/o 

CgO 



CfZ^^Zf, 



-1/2^1/2 



In order to find the /3 functions of c;,, c/, and Cg we use 
eq. ([BTT]) to find 



/3c/ = c/ 



d_ 

% 

Cb9 d 



d' 



iPf - Pi>) 



{P4> - Pb) 



fl Cge 9 



2 dg' 

JA-Pgs), (40) 

where Z = 1 + p{2/e). Finally, the /? functions read 
g^{cf - Cb) 



Gn^Cbicf + CbY 



\cf-cg){Ac) + ?,CfCg+cl 



/3(cfc) 



/3(cg 



3(4^)^Cg(Cg + C/)2 



S,TT^cjcb ' 



3(47r) 



(41) 



This is exactly expected: since only fermions can couple 
to both scalars and gauge fields, we find that the photon 
and scalars speeds of light are identical to those found 



before, while the fermion speed gets contributions from 
both Yukawa and gauge sectors. 

We can see that c/ — Cb — Cg is an IR attractive fixed 
line by computing the eigenvalues of the Jacobian J = 
{dfii/dcj) |c/=c6=cj,, where i,j = Cf,Cb,Cg, 



(7g^ + 12e2) ± J(7g^ + 12e^? ~ 30Ag^e^ 



(42) 



6.2. The general case 

We consider Nf fermions interacting with Nb scalars 
or gauge bosons. Although we shall carry out the cal- 
culation in the case of Yukawa interactions, the abelian 
and non-abelian /3-functions have the same structure as 
we pointed out above. 

The general Lagrangian density reads 

£ = iiJal^dotpa - iCf^lpal ■ d'llJa + ■^do<l)ido(j)i (43) 

where summation over repeated indices is implied. De- 
noting — Uab + ivl^b ^^'^ noticing that z' — > z't as we 
move z' across the vertex [g] we find 



Zf=l- 



Zrh. = 1- 



2 Y-«a(c/.+2c6.) /2 



3(47r)^^ Cf^cb^icf, + Cb,)"^ \ej ' 

16 'a,4: fl' 



8 (4.+4c/.c/,+c2j 

^b, = 1 "T / , — 



(44) 



Zf^, Z^^ can be read directly from the Yukawa expres- 
sions in Eq. [T71 while Z^. and Zb^ are obtained using a 
series of integrals similar to those given in Appendix A. 
Notice that quantum loops can also generate ofii'-diagonal 
corrections Z^^^^^ if the couplings contain off-diagonal 
components. These corrections will induce kinetic mix- 
ing terms of the form iuabi^ad'^j'^'ipb + iPab4'ad'''^^'4>b- In 
Lorentz invariant theories, where Uab = Pab, we can al- 
ways find basis where aab = Pab are diagonal by per- 
forming SO{Nf) rotations. However, in the present case, 
and since in general Uab l^ab, we have the freedom to 
diagonalize either the time-time or the space-space com- 
ponents. Diagonalizing the time-time component, and 
hence working in a basis where we have canonical kinetic 
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terms, will always leave space-space mixing terms. We 
assume that these terms are always small compared to 
the diagonal speeds, i.e. /3a&/c/„ << 1 for all a, and 
we ignore their evolution in the following analysis. The 
same thing can also be said about kinetic mixing terms 
for bosons. 

To be able to use the grand formula (|B11|) we define 

Z'ab = 93tab , Z*l = gg^f, , C/„ = gia , = g2t ,(45) 



and 



Zci>i — Zii , — Z2i , Z^i — Z[ 



3iab ■ 



Zsiab I 



z 



5a ■ 



Now, we write (jBllI) as 

l3f_iM = ^gf^RI ^ TlyN^^M ^ Pa09aO 



aO 



dg. 



'aO 



(46) 



(47) 



where the Greek incices run from 1 to 3 and 3, and 
the upper case lattin letters run over a and i. The non 
zero values in ([47]) are nii',2i — —Sii'/2, n2i' ,2i — /2, 

riia'M = ^aa'y '^5a^la = -Saa', P3,abi = 1/2, and 

P3,abi — 1/2- Hence, we obtain 
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FIG. 2: The smallest eigenvalue of J (vertical) against the 
number of fermions Nf (horizontal) taking A'^;, = 1. The cal- 
culations are based on values of 2^6 between and 1 generated 
randomly at each Nf. The plot shows that A oc Nf for large 
Nf. 



IMPLICATIONS OF LOGARITHMIC 
RUNNING 



PicbJ 



3 (47r)' Cb, 



■^ab^ba 



a.b 



(c/„ 



3 (4-)' ir 



C&i (c/i, + Cbi) 



(48) 



Again, we find in this general setup that c/^ = Cb^ = c 
for all i and a is a fixed line. To study the nature of this 
line we perform a perturbation to the system about this 
line, i.e. we construct the Jacobian matrix at c/^ = c;,. = 



Defining AJ^ 



and K = J: 



J, 



(Nf+Nb)x(Nf+Nb) 



-^NfXNf 
■^NbXNf 



b^a ab ba 



ztl we find 



(49) 



where 



/I 



3 (47r)' c3 
-2 [K + K] 



3 (47r)' c3 



' 3 (47r)' c3 

-4[A'. + r^] 

(4^)' c3 ' 



J, 



2J 



4 Eg K + r^] 

(47r)^ c3 



(50) 



Although we were not able to diagonalize J analytically, 
numerical calculations show that we always have a spec- 
trum of positive eigenvalues on the top of a zero mode. 
In figure ^ we plot the smallest eigenvalue A of J, which 
governs the behavior of c/^ and c^. , against the number 
of fermions NS ^ for fixed number of bosons f^b- We find 
that A oc TVy for large Nf. Hence, the effective running 
increases with the number of species as expected. 



Let us consider how the differing speeds approach each 
other by studying the situation where the speeds are rel- 
atively close, but not identical, at some scale /u*. We 
treat this problem to first order in the speed difference. 
If we define the relative speed difference as 



Cb , 

77 = 1 , 

the generic beta function has the form 



drj bg^ 



77 + 0(77'), 



(51) 



(52) 



where 6 is a constant of order unity and where we denote 
the common low energylimit as Cf ~ Cb ~ c. 

If the coupling g were to be treated as a constant, this 
running could be integrated to yield 



r/(/i) = ?7, 



(53) 



The rescaling is power-law in form. If the coupling con- 
stant is small, the power-law exponent is also small and 
the running is slow. However, if the coupling is large (and 
constant) the running would be rapid with a power-law 
form, leading quickly to a universal speed of light. 

However, the coupling itself also runs. For example, 
the Yukawa coupling beta function can be integrated to 
yield 



47rg2 



W 51og(A^) 



(54) 



with A > fi. While this coupling could be large at high 
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energy it runs to smaller values at low energy^. This 
produces a quite different form for the running of the 
relative speeds. The correct form for the running of rj is 



log 



lof 



(55) 



. This implies that the difference in the speeds runs only 
logarithmically towards each other. 

There are tight constraints on the equality of the lim- 
iting velocities for the different particles. For direct mea- 
surement of the velocities, we can look at timing accom- 
plished at high energy accelerators. For example at LEP, 
the electron beam travels at essentially the limiting ve- 
locity, since E/m = 7 ss 10^. The timing of the ac- 
celerator relies on this limiting velocity being the speed 
of light. Because the timing of each bunch is recorded 
within ±50 ns over about 1000 revolutions in the 27 km 
accelerator |S|, we estimate that this constrains 77 < 10~^ 
for electrons. 

However, indirect constraints are more powerful, and 
these have been described by Altschulfl0|. For Ce > c, 
energetic electrons traveling faster than the speed of light 
will radiate Cherenkov light, losing energy until they 
move at only the speed of light. This effect produces 
a maximum energy for subluminal motion, which is con- 
strained by the observation of energetic electrons in astro- 
physics. For Ce < c, there is a constraint from the cutoff 
frequency in synchrotron emission. These constraints are 
more powerful than direct measurements because they 
bound factors of 7c = 1 / ■\/ 1 — /c^ « rather than 

the linear bounds on rj from the velocity measurements. 
Altschul's bounds are \ri\ < 10~^^. 

In order to achieve this close equality of the different 
speeds with logarithmic running, the running needs to 
occur over an exponentially large energy range. For ex- 
ample, even if we take 77, ~ 10~^ and A/^* ~ 2 (which 
barely allows perturbation theory to be used near the 
energy /i*), we would need log(A/me) ^ 10^"^, where we 
have generously used rrie as the low energy scale. This 
clearly poses a problem for model building. 



8. EMERGENT LORENTZ SYMMETRY: 
TOWARD MODEL BUILDING 

In any realistic model of emergence without an intrinsic 
Lorentz invariance, we do not expect the different species 
to emerge with the same limiting speed. In the above 
sections we showed that there is a potential mechanism 



to overcome this problem in a class of models whenever 
we run the renormalization group down to lower energies. 
However, we found that the speeds of light are forced 
to run logarithmically along with the running coupling 
constants. This is a relatively slow running if we want 
to meet the stringent constraints on Lorentz violations 
without having to fine tune the speeds at the UV. In this 
section, we propose a way out of this situation. 

In order to increase the effect of RG running there 
are two options. One is to keep the coupling constant 
large and unchanged with energy scale. Such a nearly 
conformal theory would convert logarithmic running into 
power-law running, as we saw in the last section. We 
also need the large coupling such that the exponent is 
large. Such theories are under active investigation [Tlj in 
the context of "walking Technicolor" where slowly run- 
ning but strongly interacting theories are used to provide 
dynamical breaking of the Electroweak Theory while not 
producing excessive flavor changing processes. Should 
walking Technicolor theories prove successful, it would 
be quite interesting to tie those results with the idea of 
an emergent limiting velocity. The other option is if there 
are a very large number of fields of different scales, such 
that the running is increased by a large (and energy de- 
pendent) factor. We explore this option below. 

We introduce a large number Nf of hidden fermions in 
addition to the Standard Model (SM) ones T?]. More- 
over, we assume that all these fermions (hidden and SM) 
have the same origin, and hence all have the same ini- 
tial speed of light 1 -f- c/, , with |c/, | << 1 , at some 
UV emergence scale /i*. As a warm up calculation, we 
assume that the fermions have a common initial charge 
e* under a single U{1) gauge sector. The gauge photon 
emerges with some initial speed 1 -I- Cg^ , with jc^J << 1 
, that is different from the speed of fermions. At the UV 
scale the fermions are taken to be massless and hence 
will participate in the running of the gauge coupling as 
well as photon and fermions speeds. As we run down 
our RG equations, some of the hidden fermions become 
massive and decouple from the RG equations. We model 
the dependence on the mass scale using a power law 



(56) 



where F/ and a/ are positive constants, and Mf is an 
IR mass scale. Since the fermions have a common 
initial speed and a common initial coupling strength, the 
evolution of the system can be modeled with a single c/ 
and e common to all fermions. Therefore, to one-loop 



Clearly for the running coupling in Yang-Mills theory, the cou- 
pling is small at high energies and becomes large only at low 
energy. However the essential point is the same - that the cou- 
pling constant does not remain large at all energy scales. 



^ This exact behavior is also exhibited in models of large extra- 
dimensions where the Kaluza-Klein (KK) modes (from the 4D 
point of view) obey a power law as in Eq. I56lll3ll . In this context, 
Mf is the lowest KK mode Mf ~ 1/L, where L is the size of the 
extra dimension, and a = d is the number of extra dimensions. 
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order we have 



de 

1^ ~ 127r2' 
dcf _ e^{cf-Cg) 



dfi 

dCn 



37r2 

e^(Cg - Cf) 



Nfi^,). 



dn Gtt^ 
Integrating this system yields the running charge 



1 + 



\Mf J { Mf ) 



(57) 



(58) 



and speeds 



(59) 



Hence, we see that both e{p) and Cg{fi) — Cf{fi) experience 
power-law running with IR values given by 



,(0)-c^(0) ^ e^{0) _ QiT^af 



2r 



/ 



My 



(60) 



Therefore, we can choose fi^/M ~ lO^**/"/ in order to 
meet the stringent requirement rj ~ 10~^*. However, 
taking < 1, so that we can trust our perturbation 
theory, weakens the coupling strength to values ^ 10~^^. 
This is way below any interesting phenomenology. 

In order to cure this problem, we introduce a large 
number of hidden U{1) sectors in addition to the SM 
U{l)hyp- We also assume that all these gauge sectors 
emerge with the same initial speed of light. In addition, 
we take all fermions to be charged under the different 
C/(l)s with the same initial charge. As in the case of 
fermions, we assume that the hidden J7(l)'s are massless 
at the UV scale, then they become massive and decouple 
as we run down the RG equations. Hence, the gauge 
sector obeys the scale-dependent relation 



Under these assumptions the second Eq. 
placed by 



(61) 



in (1571) is re- 



de/ 
dfi 



37r2 



(62) 



The solution of e(/i) is still given by Eq. [58l while that 
of Cg(/i) — c/(/i) can be expressed in terms of the hy- 
pergeometric function. Instead, we numerically integrate 
our system setting appropriate values of the parameters 
af, ttg, Tf, and Tg. From FIG. [3] we see that introduc- 
ing many copies of U{l)s achieves power-law suppression 



0.1 



10- 



10- 



10- 



10- 




FIG. 3: Numerical simulation of the running of the charge 



e, and ratio 77 

of log(M/M/) 
-2 



- Cf/Cg. 

We take M„ 



The horizontal axis is in units 



Mf 



1, To 



lo-^^/ = 5, 

Of = 10"^, fJ.*/Mf = 10", and Og = 1, 1.2, and 2 for the 
dashed, dot-dashed and continuous lines respectively. We also 
use the initial conditions c/, = 0.6, Cg, = 0.3, and e, — 0.5. 
The running of e is logarithmic, while the running of rj is 
power-law. Very small values of 77 are achieved in a very 
short interval of running as we increase ag, i.e. as we increase 
the number of gauge sectors. We also note that for Og = 2 
the value of rj is saturated by the error tolerance of the code. 
More powerful computations may give smaller values. 



of ?; in a very short interval. Since the many U{l)s do 
not intrude the evolution of e, we can still get reasonable 
coupling strength in the IR. 

To understand the choice of parameters used in the 
simulation in FIG.[3l it is instructive to calculate the to- 
tal number of fermions and gauge fields as seen in the 
UV. Using Eqs. [521 and [HIl and the numerical coeffi- 
cients given in FIG. |3l (take ag = 2) we find Ng w 10^, 
and Nf « 6, which explains the above findings. Because 
the RG fiow of e is enhanced only by a few number of 
fermionic species, the coupling constant runs only loga- 
rithmically. While the huge number of gauge fields that 
participate in the RG equation of speeds force the run- 
ning oi Cg — Cf to be extremely fast. 

Because the freezing of the speeds sets in almost im- 
mediately, we can try to replace the power law numbers 
in Eqs. I5S1 and ICTl bv constant numbers. The idea is that 
we just need to have Ng » Nf » 1 for a short interval 
until the speeds freeze to their desired values. Then, we 
can immediately integrate the RG equations to find 



e'(M) 
Cg. - Cf, 



Nf 



■l°g(lf) 



(63) 



Now if we take e* « 1, ft*/ fiiR. 
Ng sa 900 we find that e]J^TT r 
Moreover, assuming that the masses of the hidden sectors 
are larger than fim, these masses decouple below /i/^ and 



192, Nf w 120, and 
1/129, and rj < 10~". 
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drop from the RG equations. Therefore, one needs only 
a constant large number of copies, Ng ,Nf » 1 to ac- 
complish the emergence of an IR Lorentz-invariant fixed 
point in a short interval of running. Moreover choos- 
ing the ratio Ng/Nf >> 1, we can meet the stringent 
constraints on the parameter rj. This opens up the pos- 
sibility that many copies of hidden sectors may suppress 
Lorentz- violating effects already present at the TeV scale. 



9. DISCUSSION 

Achieving a universal speed of light is a challenge for 
theories which do not postulate an fundamental Lorentz 
symmetry. This problem is visible in known emergence 
models Q and also in Hofava-Lifshitz gravity [Ij. For 
emergent gauge fields, the Weinberg- Witten theoremQ 
suggests that this will be a continual challenge as a 
Lorentz-noninvariant initial theory may be required. 

We have shown through several examples that a com- 
mon limiting velocity can be emergent at low energies 
even if the original high-energy theory involves fields sat- 
isfying the wave equation with different velocities. There 
is a heuristic rationale for this in that since fields can 
transform into each other through interactions, the end- 
point where all the fields travel in unison is preferred. 
The renormalization group treatment indeed produces 
this outcome. 

Because the running is only logarithmic for simple sys- 
tems, it would take an exponentially large amount of run- 
ning in order that the limiting velocities be close enough 
to agree with experiment. We addressed the phenomeno- 
logical constraints in Sec. 7. However, power law running 
is also possible if the coupling is large and constant, or 
if there are a very large number of interacting degrees of 
freedom. We have reported on a model with this latter 
property. 

It is important to note that not all forms of Lorentz- 
violation disappear at low energies. The renormaliza- 
tion of a general parameterization of Lorentz-violation 
of QED has been studied in Ref. Q, and some opera- 
tors that grow at low energy are found. A well behaved 
emergent theory must avoid those operators. 

The running of the limiting velocities only happens due 
to the interactions that couple one type of particle to an- 
other. This implies that the running will be weakest if 
the coupling is weak. At low energies the gravitational 
coupling is by far the weakest of all the fundamental 
forces. This implies that the most plausible velocity dif- 
ference would be that of gravity. While there have been 
some claims that the speed of gravity has been indirectly 
measured [14], the consensus appears to be that there is 
no experimental constraint on the speed of gravitvfisj. 
However indirectly there is a stringent limit at the lO"'^^ 
level on the difference of the speeds of gravity and that of 
light from gravitational Cherenkov radiation jlGj which 
is valid if the speed of gravity is less than that of light. 
Future experiment with gravitational wave detectors pro- 



vide the best opportunity to measure or constrain the 
difference if the speed of gravity is greater than that of 
light. 
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Appendix A: Useful Integrals 

In this appendix we work out the details of the Inte- 
grals I*^ and / appearing in Eq. |51 These integrals are 
given by 



/o = / ^fcO 
27r 



dk 



1 



{2ny (fcO)2 - fc2 
1 



{p"/cb - Cfk"/cbf - [k- p 



(Al) 



and 



dk 



k 



dk° 

^7 (27r)^ (fc0)2 - fc2 
1 

(p°/c6 - Cfk° /cb)^ - {k-p^ 



(A2) 



To perform the integral 7*^, wc first use the Feynman trick 
to find 



dx- 



1 



dk 



{2^ f Jo (fc2„A2)" 



(A3) 



where A2 = -x{l - x)p^ + x{k^f + (1 - x){p°/cb - 
Cfk^/cb)'^. Next, we interchange the integrals dx and 
dk and perform the integral over dk to find 



/0 = 



(4^) 



—rr- 

3/2' I 2 



dk 



1 



„ k° / dx- , 



(A4) 



Further, we exchange the integrals dx and dk^ , and rear- 
range the integrands to find 



/o = - 



(47r 



>3/2^ \ 2 



dx- 



1 



^:r(l~c2./c2)+c2./ 



27r 



^kl + 2fcoi?o - 



(A5) 
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where 



(1 - x)p°Cf/cl 



x{l-c)/cl)+c}lcr 
x{l~c)/cl)+c)lcl 



(A6) 



Then, we perform the integral over dk^ , after analytically 
continuing from 1 to c? = 1 — e dimensions, to obtain 



(4^) 



2 „2i 



dx- 



x{l - c)/cl) + cj/cl 

r(6/2) 



3/2- 



e/2 



x{fYlcl - x f x + (1 - a;)c2/c2 



Finally, we find 



2C6 



finite 



(47r)^c^(l + a)^ Ve 
where a = Ch/cf. Similarly, we can show 

/=^Ml±Mf2^fi,ite) 
(47r)' 3(l + a)2 \t J 



e/2 



(A7) 



(A8) 



(A9) 



The vertex correction in Eq. [Tl] results in the integral 



igG 



dK {K' 



0\2 



{2TTf c){K°Y/cl-K^ 




(AlO) 



Next, we use The Feynman trick to get 
2g3 f dK" f dK 



-igG = 



27r 



1 nl-x 

dx / dy— 
Jo 



3 yK^-iK°y 



1 



(All) 



where 



if2_ A2 



[l + y(-l + 4/c?)] (K^ 
-2[{l-x-y)P?+P°x] K° 

+ (1 - X - y)Pi + xP2 
+ {l-x~y)Pl+xPi. 



(A12) 



Then, proceeding as we did before, we finally obtain the 
result in EalT^ 



Appendix B: A general Setup to calculate the /3 
functions 



We assume that the parameter space is spanned by 
gi, i,j = 1,2,..C couplings (these could be coupling 
strengths as well as speeds). Quantum loops will gen- 
erate Zm, m,l = I, 2, ...D corrections to the original La- 
grangian, and we restrict our treatment to one loop order. 
In general, we may write 

D 

9^0 = g^M U ZZ-Hf^)f^''" ■ (Bl) 
m— 1 

Taking the derivative of of (jBip w.r.t. fi we obtain 

D 



m—l 
D 



D 



1^1 



(B2) 



Writing Z„,{p) = 1 + p.mi^^)^ we find Z,'(^) = 
Also, using the definition A(/x) = 

eq. (|B2|) becomes 



m—l 

D C 



1 = 1 j=l ^ m^l 



(B3) 



Since we are only interested in one-loop corrections, we 
can ignore all terms O (l/e^). Hence, eq. (jB3l) reads 



1 + - E ^™(^)^ 



m—l 



D 



1 + - ^ Pm(M)»' 



m—l 



. (B4) 



Now, we can rewrite eq. (IB4p in the following simple 
expression 



l + -A^{^i) 



-^C,y(Ai)/3j(^) 



D 



(B5) 
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where 



Hence, solving for /3 from (|B7p we obtain 



D 



dgi 



D 



(B6) 



m— 1 

Eq. (IBSp can also be written in the matrix form 

D 



M (a^) P {p) = - 9P 



m— 1 



(B7) 



where 



X (m) = 



iC21 1 + 1^2 

The inverse of M. (m) is given by 



7C2C 



.(B8) 



« -1 
M (m) 



/1-|A -|Cl2 

IC21 1-M2 

+0- ' 



(B9) 



\ m=l J 



2EC,j3,Pj.(B10) 



Finally, we rearrange the terms to find 



ft = E E 



(Bll) 
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